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Abstract 

We present an improved measurement of CP-violating asymmetries in B° — > tt + tt~ decays 
based on a 78 fb -1 data sample collected at the T(4S) resonance with the Belle detector at the 
KEKB asymmetric-energy e + e~ collider. We reconstruct one neutral B meson as a B° — > tt + tt~ 
CP eigenstate and identify the flavor of the accompanying B meson from inclusive properties 
of its decay products. We apply an unbinned maximum likelihood fit to the distribution of the 
time intervals between the two B meson decay points. The fit yields the CP-violating asymmetry 
amplitudes Attt = +0.77 ±0.27(stat) ±0.08(syst) and S nn = -1.23 ±0.41(stat) Ig^syst), where 
the statistical uncertainties are determined from Monte Carlo pseudo-experiments. We obtain 
confidence intervals for CP-violating asymmetry parameters Ann and S n n based on a frequentist 
approach. We rule out the CP-conserving case, Ann = $nn = 0, at the 99.93% confidence level. 
We discuss how these results constrain the value of the CKM angle 4>2- 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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I. INTRODUCTION 



In 1973, Kobayashi and Maskawa (KM) proposed a model where CP violation is accom- 
modated as an irreducible complex phase in the weak-interaction quark mixing matrix [ll. 
Recent measurements of the CP- violating parameter sin 2(f) 1 by the Belle and BaBar ^| 
collaborations established CP violation in the neutral B meson system that is consistent 
with the KM model. Measurements of other CP-violating parameters provide important 
tests of the KM model. 

The KM model predicts CP- violating asymmetries in the time-dependent rates for B° 
and B° decays to a common CP eigenstate, fcp 0]. In the decay chain T(4S') — > B°B° — > 
fcpftug, in which one of the B mesons decays at time tcp to fcp and the other decays at 
time t tag to a final state that distinguishes between B° and B°, the decay rate has a 
time dependence given by p 

e -|At|/r B o 

V« n (At) = — [1 + q- {S nw sin(Am d At) 

At b o 

+A wn cos(Am d At)}], (1) 

where t#o is the B° lifetime, Arrid is the mass difference between the two B° mass eigenstates, 
At = tcp — ttag, and the 6-flavor charge q = +1 (—1) when the tagging B meson is a B° 
(B°). The CP- violating parameters «S 7r7r and A n7T defined in Eq. ((TJ) are expressed as 

|A| 2 - 1 2Im\ 

™ ~ |ApTT' ~ |ApTT' 

where A is a complex parameter that depends on both B°-B° mixing and on the amplitudes 
for B° and B° decay to 7r + 7r~. In the Standard Model, to a good approximation, |A| is 
equal to the absolute value of the ratio of the B° to B° decay amplitudes. A measurement 
of time-dependent CP- violating asymmetries in the mode B° — > tt + tt~ p] is sensitive to 
direct CP violation and the CKM angle 02 @- If the decay proceeded only via a b — > u 
tree amplitude, we would have S n7T = sin20 2 and A^n = 0, or equivalently |A| = 1. The 
situation is complicated by the possibility of significant contributions from gluonic b —>■ d 
penguin amplitudes that have a different weak phase and additional strong phases As a 
result, S nn may not be equal to sin2</> 2 and direct CP violation, A n7T ^ 0, may occur. 

Belle's earlier published study [l(j was based on a 41.8 fb _1 data sample containing 44.8 
x 10 6 BB pairs produced at the T(4S') resonance. The result suggested large direct CP 
asymmetry and/or mixing- induced asymmetry in B° — > 7r + 7r~ decay while the correspon ding 
result by the BaBar collaboration based on a sample of 88 x 10 6 BB pairs did not [11]. 
In this paper, we report an updated measurement that is based on a 78 fb _1 data sample, 
containing 85 x 10 6 BB pairs. The most important change is that we determine the statistical 
significance and uncertainties in the CP parameters from the distributions of the results 
of fits to Monte Carlo (MC) pseudo-experiments, instead of from errors obtained by the 
likelihood fit to experimental data. In addition, we have made three significant improvements 
to the analysis: a new track reconstruction algorithm that provides better performance; a 
new proper-time interval resolution function that reduces systematic uncertainties; and the 
inclusion of additional signal candidates by optimizing the cuts for continuum background 
suppression. 

In Sec. |H] we describe the KEKB collider and the Belle detector. The reconstruction 
of B° — > 7r + 7r~ decays is described in Sec. IIHI The flavor-tagging procedure and vertex 
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reconstruction are described in Sees. |lV|and|V| After discussing the signal yield in Sec. EU 
and introducing the method to measure A nn and S nn from At distributions in Sec. I VI II 
we present the results of the fit in Sec. I VI III and discuss constraints on 2 in Sec. IIXI We 
summarize the results in Sec. OD 



II. EXPERIMENTAL APPARATUS 

The data reported here were collected with the Belle detector at the KEKB asymmetric- 
energy e + e~ collider 0], which collides 8.0 GeV e~ and 3.5 GeV e + beams at a small (±11 
mrad) crossing angle. The T(4S) is produced with a Lorentz boost of (3^ = 0.425 nearly 
along the electron beamline (z). Since the B° and B° mesons are approximately at rest in 
the T(4S) center-of-mass system (cms), At can be determined from Az, the displacement 
in z between the fcp and / tag decay vertices: At ~ (zqp — Zt^/P^c = Az/P^c. 

The Belle detector is a large-solid-angle general purpose spectrometer that consists of a 
silicon vertex detector (SVD), a central drift chamber (CDC), an array of aerogel thresh- 
old Cerenkov counters (ACC), time-of- flight scintillation counters, and an electromagnetic 
calorimeter comprised of CsI(Tl) crystals located inside a superconducting solenoid coil that 
provides a 1.5 T magnetic field. An iron flux return located outside of the coil is instrumented 
to detect K Q L mesons and muons. For more details, see Ref. \l~ 



III. RECONSTRUCTION OF B° -> vr+vr DECAYS 



The B° — > 7r + 7r event selection is described in detail elsewhere 13 . We use oppositely 



charged track pairs that are positively identified as pions according to the likelihood ratio for 
a particle to be a K ± meson, KID = £(K) /[£(K)+£(ir)], which is based on the combined 
information from the ACC and the CDC dE/dx measurements. Here we use KID<0.4 as 
the default requirement for the selection of pions. For tracks in the momentum range that 
covers the B° — > tt + tt~ signal, this requirement has a pion efficiency of 91% and 10.3% of 
kaons are misidentified as pions ( 10.0±0.2% from K~ and 10.6±0.2% from K + ). 

We also select B° — > K + tt~ candidates, which have the same track topology as B° — > 
7r + 7r~ candidates, by positively identifying one charged track as a kaon and the other as a 
pion. We use KID>0.6 for the selection of kaons. This requirement has a kaon efficiency of 
84% and a misidentification rate from pions of 5%. 

Candidate B mesons are reconstructed using the energy difference AE = E^ 113 — E^ m 



and the beam-energy constrained mass M bc = y (E^ m ) 2 — (p^) 2 , where E™* m is the cms 
beam energy, and E^ and p^° s are the cms energy and momentum of the B candidate. The 
signal region is defined as 5.271 GeV/c 2 < M bc < 5.287 GeV/c 2 and \AE\ < 0.057 GeV, 
corresponding to ±3a from the central values. In order to suppress background from the 
e + e~ — > qq continuum (q = u, d, s, c), we form signal and background likelihood functions, 
£s and £bg, from two variables. One is a Fisher discriminant determined from six modified 
Fox- Wolfram moments |l5j ; the other is the cms B flight direction with respect to the z axis. 
We determine £$ from a GEANT-based Monte Carlo (MC) simulation [16(, and £bg from 
sideband data in the 5.20 GeV/c 2 < M bc < 5.26 GeV/c 2 and -0.3 GeV < AE < 0.5 GeV 
region. We reduce the continuum background by imposing requirements on the likelihood 
ratio LR = £ s /(£ s + £bg) f° r candidate events, as described below. 
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IV. FLAVOR TAGGING 



Leptons, kaons, and charged pions that are not associated with the reconstructed B° — ► 
7T + 7r _ decay are used to identify the flavor of the accompanying B meson. We apply the 
same method used for the Belle sin 20! measurement 0]. We use two parameters, q and r, 
to represent the tagging information. The first, q, is defined in Eq. (jTJ). The parameter r is 
an event-by-event, MC-determined flavor-tagging dilution factor that ranges from r = for 
no flavor discrimination to r = 1 for unambiguous flavor assignment. It is used only to sort 
data into six r intervals. The wrong tag fractions for the six r intervals, Wi (I = 1,6), are 
determined from the data and are summarized in Table UJ 

TABLE I: The wrong tag fraction wi for each r interval. The errors include both statistical and 
systematic uncertainties. 



I r interval wi 

1 0.000 - 0.250 0.458±0.007 

2 0.250 - 0.500 0.336±0.010 

3 0.500 - 0.625 0.228±0.011 

4 0.625 - 0.750 0.160±0.014 

5 0.750 - 0.875 0.112±0.015 

6 0.875 - 1.000 0.020±0.007 



In the previous publication [lOj . we required LR > 0.825 for all candidate events, while in 
this analysis we optimize the expected sensitivity by including additional candidate events 
with a lower signal likelihood ratio. The requirements on LR vary for different values 
of tagging dilution factor r, as indicated in Table |H1 since the separation of continuum 
background from the B signal varies with r; there are 12 distinct regions in the LR-r plane. 

V. VERTEX RECONSTRUCTION 

The vertex reconstruction algorithm is the same as that used for the sin20! analysis 0]. 
The vertex positions for the fcp decay (7r + 7r~) and / tag decays are reconstructed using tracks 
with associated hits in the SVD. Each vertex position is also constrained by the interaction 
point profile, determined run-by-run, smeared in the r-<p plane to account for the B meson 
decay length. With these requirements, we are able to determine a vertex even with a single 
track; the fraction of single-track vertices is about 10% for zqp and 22% for z taug . The / tag 
vertex is determined from all well-reconstructed tracks, excluding the B° — > 7r + 7r~ tracks 
and tracks that form a Kg candidate. 

VI. SIGNAL YIELD 

Figures^a) and (b) show the AE distributions for the B° — > 7r + 7r~ candidates that are in 
the Mb c signal region with LR > 0.825 and with LR < 0.825, respectively, after flavor tagging 
and vertex reconstruction. In the Mb c and AE signal region, we find 275 candidates for LR 
> 0.825 and 485 candidates for LR < 0.825. The B° -> signal yield for LR > 0.825 is 

extracted by fitting the AE distribution with a Gaussian signal function plus contributions 
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AE (GeV) AE (GeV) 



FIG. 1: AE distributions in the Mb c signal region for (a) B° — > tt + tt candidates with LR > 0.825, 
(b) B° -> tt+tt- candidates with LR < 0.825, (c) L?° -> candidates with Li? > 0.825, and 

(d) 5° -» K+vr- candidates with Li? < 0.825. The sum of the signal and background functions 
is shown as a solid curve. The solid curve with hatched area represents the tt + tt~ component, 
the dashed curve represents the K + tt~ component, the dotted curve represents the continuum 
background, and the dot-dashed curve represents the charmless three-body B decay background 
component. 

from misidentified B° — ► K + tt~ events, three-body 5-decays, and continuum background. 
The fit yields 106±i| tt+tt events, 4llg° K+tt events and 128^g continuum events in the 
signal region. The errors do not include systematic uncertainties unless otherwise stated. 
Here the error on the yield of continuum events in the signal region is obtained by scaling 
the error of the yield from the fit that encompasses the entire AE range. For LR < 0.825, 
we fix the level of tt+tt^ signal by scaling the LR > 0.825 number by a MC-determined 
factor and that of the continuum background from the sideband. The ratio of the K+tt^ 
background to the tt + tt~ signal is fixed to the value measured with the LR > 0.825 sample. 
We obtain 57 ± 8 tt + tt~ events, 22lJ K + tt~ events and 406 ± 17 continuum events in the 
signal region for LR < 0.825. The contribution from three-body -B-decays is negligibly 
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small in the signal region. Figures [IJc) and (d) show the AE distributions for the selected 
B° — > K + 7r~ candidates. 

VII. MAXIMUM LIKELIHOOD FIT 

The proper-time interval resolution function R wn for B° — > 7r + 7r~ signal events is formed 
by convolving four components: the detector resolutions for zqp and z tag , the shift in the 
-Stag vertex position due to secondary tracks originating from charmed particle decays, and 
the smearing due to the kinematic approximation used to convert Az to At. We use the 
same parameters as those used for the sin20! measurement We determine resolution 
parameters from fitting the data for the neutral and charged B meson lifetimes. A small 
component of broad outliers in the Az distribution, caused by mis-reconstruction, is rep- 
resented by a Gaussian function. The width of the outlier component is determined to be 
42^4 ps; the fractions f \ of the outlier components are (2 ± 1) x 10~ 4 for events with both 
vertices reconstructed with more than one track, and (2.7 ± 0.2) x 10~ 2 for events with one 
or two single-track vertices. We assume R nn = Rk-k and denote them collectively as R s i g . 
The parameters of the continuum background resolution function R qq are obtained from the 
sideband data. 

The CP asymmetry parameters, A nn and S niT , are obtained from an unbinned maximum 
likelihood fit to the observed proper-time distribution. For this purpose, we use probability 
density functons (PDFs) that are based on theoretical distributions that are diluted and 
smeared by the detector response. The PDF for B — > 7r + 7r~ signal events (V^n) ls gi ven by 
Eq.((TJ), with q replaced by q(l — 2wi) to account for the dilution due to wrong flavor tagging. 
The PDF for B° -> K+n- background events is V q K7T {At, w t ) = e^'l^ /(4r B o){l +q-{l- 
2wi) ■ Ak-k ■ cos (Am d At)} , where we assume as a default that there is no CP asymmetry 
for the B° — > K + tt~ mode. The effect of a possible non-zero value for Ak-k is determined 
by varying Ak-k by the error obtained from fits to the self-tagged B° — > K + tt~ sample 
and is included in the systematic error. The PDF for continuum background events is 



V q q(At) = (1 + g-Ak g )/2{/ r e-l A 'l/^/(2r b kg) + (1 - f T )6{At)}, where f T is the fraction of 



the background with effective lifetime r^ g and 5 is the Dirac delta function. For B° — > 
7T + 7r _ candidates where both vertices have at least two tracks, we use f T = 0.014 lo!oo4 an d 
2.37« J! ps, which are determined from the events in the gg-background-dominated 
sideband region: 5.20 GeV/c 2 < M bc < 5.26 GeV/c 2 and 0.10 GeV < AE < 0.50 GeV. 
For events with a single-track vertex, we use f T = 0. The effect of the uncertainty in ^4bkg, 
determined by varying Ahkg by the error from the fit to the sideband data, is included in 
the systematic error. 

We define the likelihood value for each (i-th) event as a function of A^ and S W7T : 



Here the probability functions f™ (k = nir, Kir or qq ) are determined on an event-by- 
event basis as functions of AE and Mb c for each LR-r interval (m— 1, 12). For example, 



fZ(AE,M bc ) is F^J(F^ + F W7r g^ + F K7r g^), where g£ is the average fraction of 



event-type k for the m-th LR-r interval (g™^ + g^ w + g"± = 1). We determine these 




+f^V Kw (At', Wl )) ■ R mg {AU - At') 
+f^V q - q (At') ■ R q - q (At t - At')}dAt' + UVotiAU). 



(3) 
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parameters from the numbers of events in the sideband data and from fractions of B° — > 
7r + 7r~ MC events. Table |H] lists the values of g™ for the 12 LR-r regions. We obtain 
9Kn = 0.382 x g^ from the fit to the AE distribution for the B° — > tt + tt~ candidates with 
LR > 0.85. The distributions of AE and M^ c for the B° — > tt + tt _ signal shape function 
F nn (AE, M bc ) and B° — > K + n~ background shape function F Kw (AE,M hc ) are fit with 
Gaussian functions. F^(AE, M^ c ) is the continuum background shape function, and the 
distributions of AE and M^ c are fit with m-dependent linear functions and the ARGUS 
background function ^t}i respectively. The small number of signal and background events 
that have large At are accommodated by the outlier PDF, V a i, with fractional area / j. 

In the fit, 5^ and A^ are free parameters determined by maximizing the likelihood 
function C — fit Ph where the product is over all B° tt + tt~ candidates. 



TABLE II: The fractions of expected B — * tt + tt and continuum events for the 12 LR-r regions. 



m 


r interval 


LR interval 


m 

J 7T7T 


yqq 


1 


0.000 - 0.250 


0.825 - 1.000 


0.296 ± 0.077 


0.591 ± 0.028 


2 


0.250 - 0.500 


0.825 - 1.000 


0.385 ± 0.094 


0.468 ± 0.026 


3 


0.500 - 0.625 


0.825 - 1.000 


0.407 ±0.134 


0.438 ± 0.027 


4 


0.625 - 0.750 


0.825 - 1.000 


0.442 ±0.110 


0.389 ± 0.024 


5 


0.750 - 0.875 


0.825 - 1.000 


0.522 ± 0.081 


0.279 ± 0.022 


6 


0.875 - 1.000 


0.825 - 1.000 


0.670 ±0.129 


0.074 ± 0.009 


7 


0.000 - 0.250 


0.525 - 0.825 


0.087 ± 0.034 


0.880 ± 0.040 


8 


0.250 - 0.500 


0.525 - 0.825 


0.127 ±0.049 


0.824 ± 0.040 


9 


0.500 - 0.625 


0.425 - 0.825 


0.124 ±0.036 


0.829 ± 0.041 


10 


0.625 - 0.750 


0.425 - 0.825 


0.129 ±0.050 


0.822 ± 0.040 


11 


0.750 - 0.875 


0.425 - 0.825 


0.170 ±0.060 


0.765 ± 0.040 


12 


0.875 - 1.000 


0.325 - 0.825 


0.390 ± 0.098 


0.461 ± 0.032 



We check the validity of our fitting procedure with a large ensemble of MC pseudo- 
experiments wherein events are generated with nominal PDFs and the observed number 
of events. The parameters in the PDFs are taken from data. For various input values of 
S n7T and A n7T , we confirm that there is no bias in the fit. The MC pseudo-experiments are 
described in detail in Appendix 1X1 

VIII. FIT RESULTS 

The unbinned maximum likelihood fit to the 760 B° — > 7r + 7r~ candidates (391 B°- and 
369 £>°-tags), containing 163^23 signal events, yields A^n = ±0.77 and «S 7r7r = —1.23. 

In Figs. 12(a) and (b), we show the raw, unweighted At distributions for the 148 B°- and 127 
£>°-tagged events with LR > 0.825. The fit curves use A^ and S nn values that are obtained 
from all of the LR-r regions. The background-subtracted At distributions are shown in 
Fig. [2(c). Figure 12(d) shows the background-subtracted CP asymmetry between the B°- 
and _B°-tagged events as a function of At. The result of the fit is superimposed and is shown 
by the solid curve. 

We test the goodness-of-fit from a x 2 comparison of the results of the unbinned fit and the 
At projections for B° -> tt+tt- candidates Q. We obtain x 2 = 10.9/12 DOF (13.3/12 DOF) 
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FIG. 2: The raw, unweighted At distributions for the 275 B — > ir + ir~ candidates with LR > 
0.825 in the signal region: (a) 148 candidates with q = +1, i.e. the tag side is identified as B°; (b) 
127 candidates with q = —1; (c) B° — > TT + ir~ yields after background subtraction. The errors are 
statistical only and do not include the error on the background subtraction; (d) the CP asymmetry 
for B° — > 7r + 7r~ after background subtraction. In Figs, (a) through (c), the curves show the results 
of the unbinned maximum likelihood fit to the At distributions of the 760 B° — > 7r + 7r~ candidates. 
In Fig. (d), the solid curve shows the resultant CP asymmetry, while the dashed (dotted) curve is 
the contribution from the cosine (sine) term. 

for the At distribution of the B° (B°) tags. 

As shown in Table HTH an ensemble of MC pseudo-experiments indicates a 16.6% prob- 
ability to measure CP violation at or above the one we observe when the input values are 
A n7r = +0.57 and S wn = —0.82, which correspond to the values at the point of maximum 
likelihood in the physically allowed region (S% n + A^ n < 1); in this measurement it is lo- 
cated at the physical boundary (A^ v + = 1). This set of MC pseudo-experiments also 
indicates that for an input value on the physical boundary, the probability of obtaining a 
result outside the physically allowed region is large (60.1%). 
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TABLE III: The fractions of MC pseudo-experiments outside the physical boundary and above the 
CP violation we observe for various input values. p nir = A\ v + S% n . The selected points are on 
the line segment between {A^, S nn ) = (0,0) and (+0.57,-0.82). 



Input p nn 


The fractions outside 
the physical boundary 

(%) 


The fractions above 
the CP violation 
we observe (%) 


0.00 


1.8 


0.07 


0.20 


3.3 


0.17 


0.40 


7.3 


0.62 


0.60 


16.4 


1.7 


0.80 


34.4 


6.0 


1.00 


60.1 


16.6 



A. Statistical errors 

As described below in Section fVIIIDl we obtain confidence intervals for A n7T and S n7T with 
a frequentist approach where we use MC pseudo-experiments to determine acceptance re- 
gions, and we quote the rms values of the MC A^n and distributions as the statistical er- 
rors of our measurement. We obtain A n7T = +0.77±0.27(stat) and S n7T = — 1.23±0.41(stat). 
Here we choose values at the point of maximum likelihood in the physically allowed region, 
(A n7T , S n7r ) = (+0.57, —0.82), for the input to the MC pseudo-experiments used to obtain 
the statistical errors. The rms values determined with input values of (A^, S n7r ) = (0, 0) 
are slightly different; for these input values we obtain ±0.28 and ±0.39 for the A nn and S nn 
errors, respectively. 

In the literature, the statistical error is usually determined from the parameter depen- 
dence of the log-likelihood ratio — 21n(£/£ max ) that is obtained from the fit. Here we call this 
estimator the MINOS error, which corresponds to the deviation from the best fit parameter 
when — 21n(£/£ max ) is changed by one. The MINOS error is a convenient approximation for 
defining a 68.3% (lcr) confidence interval; however, care is needed when defining intervals 
at higher confidence levels. Figure E2 shows the log-likelihood ratio curves from our data, 
where deviations from parabolic behavior are evident; for example, 3a from the MINOS 
error for the positive side of S nn is considerably smaller than a three-standard-deviation 
error defined by the deviation from the best fit parameter when — 21n(£/£ max ) is changed 
by 9. Thus, the MINOS error overestimates the significance of a non-zero value. The 
MINOS errors obtained from the curves are also smaller than the expectations from the MC 
pseudo-experiments, as shown in Fig.0J the probability of obtaining a MINOS error smaller 
than that in our measurement is 1.2% (12.0%) for S nn {A^) [19j. These characteristics are 
reproduced in a fraction of the MC pseudo-experiments that have A n7T and S n7T input values 
that are close to the physical boundary. We therefore conclude that the rms values of the 
distributions of fit outputs of A nn and S nn from the MC pseudo-experiments, rather than the 
MINOS errors, are more appropriate as the standard statistical errors for this measurement. 
We describe an investigation of the source of the small MINOS errors in Appendix [Bj 
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FIG. 4: The result of MC pseudo-experiments with input values of = +0.57 and S n7T = —0.82: 
the distributions of (a) the negative and (b) positive MINOS errors of A n7T , and (c) the negative 
and (d) positive MINOS errors of S n7T . The arrows indicate the MINOS errors obtained from the 
fit to data. 
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B. Systematic errors 



The sources of the systematic error are listed in Table IIV1 We add each contribution in 
quadrature for the total systematic errors. We obtain 

Xtt = +0.77 ±0.27(stat) ± 0.08(syst), 
= -1.23 ±0.41(stat) +°$(syst). 

The systematic error on Ann is primarily due to uncertainties in the background fractions 
and the vertexing. For S„n, the background fractions and a possible fit bias near the physical 
boundary are the two leading components. Below we explain each item in order. 



TABLE IV: Systematic errors for Ann and <S 5 



Source 


A^n 

+error 


—error 


+error 


—error 


Background fractions 


+0.058 


-0.048 


+0.044 


-0.055 


Vertexing 


+0.044 


-0.054 


+0.037 


-0.012 


Fit bias 


+0.016 


-0.021 


+0.052 


-0.020 


Wrong tag fraction 


+0.026 


-0.021 


+0.015 


-0.016 


Physics (t b o, Am d , A K n) 


+0.021 


-0.014 


+0.022 


-0.022 


Resolution function 


+0.019 


-0.020 


+0.010 


-0.013 


Background shape 


+0.003 


-0.015 


+0.007 


-0.002 


Total 


+0.084 


-0.083 


+0.083 


-0.067 



1. Background fractions 

We estimate the systematic errors that arise from uncertainties in the parameters used 
for the event-by-event background fractions and as well as the signal fraction f^. 
Parameters that are determined from data are varied by their errors and fits are repeated; 
we add the contribution from each variation in quadrature. 

As explained in Sec. IVH[ we rely on a MC B° — > 7r + 7r - sample to determine g™, the 
background fraction in each LR-r region m. We measure the regional event fractions in 
B° — > D^tc control samples, and compare the results with those in the MC B° — > it + it~ 
sample. Each g™ value is then modified by an amount determined from the difference 
between data and MC, and from the statistical error in the control samples. We repeat the 
fit to obtain Ann and S nn , and add each difference from the nominal value in quadrature. 

The Kn background yield is obtained from the fit to the AE distribution. We estimate 
the systematic error associated with this method from an independent yield measurement 
based on a Kir enriched-sample and the Kj 7r separation performance, which will be described 

inSecEnna 

The PDF for continuum background used in the fit assumes no asymmetry (^4bk g = 0) 
between the number of events with q = +1 and with q = —1. We estimate the systematic 
error due to this assumption by varing ^4bk g by ±0.02, based on the measurement ^4bk g = 
0.013 ± 0.006 from the sideband data. 
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2. Vertexing 

We search for possible biases that may arise from the track and vertex selection by 
repeating the analysis with modified selection criteria. We include the observed changes in 
the systematic error. We also repeat the analysis by introducing charge- dependent shifts in 
the z direction artificially, and include the resulting change in the systematic error. Here the 
amount of the shift is determined from studies with cosmic rays and with the two-photon 
process e + e~ — > 7r + 7r~7r + 7r~. The systematic error associated with the IP profile is estimated 
by varying the IP smearing that is used to account for the B flight length. 



3. Fit bias and other sources 

We use large-statistics MC pseudo-experiments to determine the systematic error due 
to possible fit biases for the input A nn and S nn values near the physical boundary. We 
also perform a fit to MC B° — > 7r + 7r~ events that are generated by using a GEANT-based 
simulation. We obtain results that are consistent with input values within the statistical 
errors, which are conservatively included in the systematic error. 

Systematic errors due to uncertainties in the wrong tag fractions are estimated by varying 
each wrong tag fraction in each r region, and repeating the fit procedure. We also repeat 
the fit using wrong tag fractions obtained for B°- and 5°-tagged control samples separately. 
We add each contribution in quadrature. 

We estimate the systematic errors associated with parameters in the resolution functions, 
in the background PDF, and the physics parameters (r B o, Arrid, and Ak-k) by repeating the 
fit varying these parameters by their errors. 



C. Crosschecks 

We perform a number of crosschecks. We measure the B meson lifetime using the same 
vertex reconstruction method. The results of the application of the same analysis to various 
subsamples are also examined. In addition, we check for biases in the analysis using samples 
of non-CP eigenstates, B° — > K + n~ decays, and sideband data. 

We perform a B° lifetime measurement with the B° — > 7r + 7r~ candidate events that uses 
the same background fractions, vertex reconstruction methods, and resolution functions that 
are used for the CP fit. Figure EJa) shows the fit result. The fit to the events in the 7r + 7r~ 
sideband is also shown in Fig. EJc); the fit curve is used for the PDF of the continuum 



background. The result, t b o = 1.42lo; 12 ps, is consistent with the world-average value |21 . 

We repeat the fits for A n7T and with 7r + 7r~ candidate samples selected with more strin- 
gent selection criteria. The K + tc~ background level is reduced by tightening the accepted 
AE range or by applying more restrictive KID requirements; the continuum background is 
reduced by tighter requirements on LR and r. The effect of the At tail is checked by tight- 
ening the At range. We do not observe any systematic variation in the fit results when the 
AE, KID, LR, r, and At requirements are changed, as shown in Table IV! To account for a 
possible AE tail, we repeat the fit with an additional Gaussian function in the AE shape of 
the 7r + 7r~ signal and the K + tt~ background . The fit yields A n7T = +0.75 and S n7T = —1.21, 
consistent with our main results. In addition, we divide the data into the 42 fb _1 sample 
used for our previous measurement and the recently added sample of 36 fb -1 . The result of 



14 




At (ps) 



FIG. 5: Results of the lifetime fits for (a) 7r + 7r _ candidates and (b) K + tt~ candidates. The 
solid curves are the results of the fits, the shaded areas are the signal and the dashed curves are 
background contributions, (c) The fit to events in the 7r + 7r~ sideband. 

the new analysis on the first 42 fb _1 sample is consistent with the published result [13], and 
with that for the more recent 36 fb" 1 sample. 

A comparison of the event yields and At distributions for B°- and -B°-tagged events in 
the sideband region reveals no significant asymmetry as shown in Fig. Efa). We also use 
samples of non-CP eigenstate B° — > D~n + , D*~n + and D*~ p + decays, selected with the 
same event-shape criteria, to check for biases in the analysis. The combined fit to this 
control sample of 15321 events yields A = -0.015 ± 0.022 and S = 0.045 ± 0.033. The At 
distribution for this sample is shown in Fig. EJb). As expected, neither mixing-induced nor 
direct CP- violating asymmetry is observed. 

We select B° — > K + it~ candidates by positively identifying the charged kaons. A fit to 
the 1371 candidates (610 signal events) yields Ak-k = —0.03 ± 0.11, in agreement with the 
counting analysis mentioned above |20(, and Sk-k = 0.08 ± 0.16, which is consistent with 
zero as shown in Fig. Efc). The MINOS errors for Ak-k and Sk-k are consistent with those 
from MC pseudo-experiment models of the B° — > K + tt~ measurement as shown in Fig. [7| 
With the K + ir~ event sample, we use the vertex reconstruction method and wrong-tag 
fractions described in Sees. II VI and IVI and determine r B o = 1.46 ± 0.08 ps [Fig. Efb)] and 
Arrid = 0.55t[};Q7 ps" 1 , which are in agreement with the world average values [21(. 

The selected K + tt~ sample and the kaon mis-identification probability measured from 
a sample of inclusive D* + — > D°(^ K^-k + )-k + and <p — > K + K~ decays are used to make 
independent estimates of the K + n~ background fractions in the 7r + 7r _ sample. The results 
are 32 ± 2 K + ir~ events in the signal region with LR > 0.825 and 15 ± 2 K + tt~ events with 
LR < 0.825; these values are consistent with the results of the fit used to determine S nrr and 
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TABLE V: Selection-requirement dependence of A n7T and (MINOS errors only). 



Cut value 


A ^ 

«^*-7T7r '-''KIT 


default 

fT^TFt ^ n A \ 
^rvllJ <^ U.4J 


n 77+0.20 i oq+0.24 
u -''-0.23 1 - z,3 -0.15 


\AE\ < 2a 


n 81 +0.20 i 91 +0.25 
U - Oi -0.22 1 - Zi -0.16 
n 89 +0.21 i in+0.29 
U -° z -0.25 - i - i °-0.19 


KID < 0.20 

rvllJ <. U.1D 


n 7/ |+0.20 i n+0.26 
u -'^-0.23 1 - li -0.17 
n cn+0.22 i i 4 +0.23 

u.5y_ 24 -1.14_ 014 


LR > 0.825 
LR > 0.925 


n o/f+0.22 i iq+0.27 
U -° 4 -0.25 -l-iy_0.18 

0.69+g§g -l-24±g-fg 


\qr\ > 0.75 
\qr\ > 0.875 


1.021^ -1.2418:" 
0-9i±C;It -1.18j£* 


\At\ < 15 ps 
At < 5 ps 


n 77+O.2O i oc+0.24 
u -''-0.23 1 - zo -0.15 
n 7fi +0.20 i 9 7+0.26 
U-<D_ 22 ^'-0.17 


Sample I (42 fbr 1 ) 
Sample II (36 fb" 1 ) 


i nf) +0.l9 _i i 4 +0.30 
J..UU_ 2 5 1.14_ 021 

n Q7+0.32 _j qq+0.70 
u -°'-0.33 i - yy -0.65 
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FIG. 6: The distributions of the raw At asymmetries for (a) B° — > ir + ir~ sideband events, (b) 
the B° — > D~ir + , D*~ir + and D*~p + candidates combined and (c) B° — > K + ir~ candidates. Fit 
curves are also shown. 



A^. The changes in A^ (to^ 05 ) and S W7T (lolra) wnen these K + n~ background fractions 
are used are included in the systematic error associated with the background fraction. The 
effect of a possible charge asymmetry in the kaon misidentification rate, described in Sec. IIIII 
is negligibly small. 
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FIG. 7: The result of MC pseudo-experiments for B — > K + ir~ with input values of Ak-k = —0.03 
and Sk-k = 0.08: the distributions of (a) the negative and (b) positive MINOS errors of Ak-k-, and 
(c) the negative and (d) positive MINOS errors of Sk-k- The arrows indicate the MINOS errors 
obtained from the fit to data. 



We check the measurement of A^-n using time-independent fits to the AE distributions 
for the B° and B tags. We determine the yields from fits to the AE distributions for 
each of the 12 LR-r bins for the B° and B tags separately (24 fits in total). We obtain 
Air-* = 0. 56^^27, which is consistent with the time-dependent CP fit result. 

As discussed above, the nominal fit result is outside of the physical region. We also 
consider fits that constrain the results to be in the physical region denned by A^ n + S^ n < 1. 
The disadvantage of the constrained fitting method is that when the fit result is close to 
the physical boundary, the errors returned from the fit are not Gaussian and are difficult 
to interpret. A constrained fit finds A^ = +0.57 and «S 7r7r = —0.82, on the boundary of 
the physical region; \ 2 values that are defined in Sec. IVIIII for the At projections are x 2 = 
12.4/12 DOF (13.6/12 DOF) for the B° (B°) tag. 



D. Significance 

We use the Feldman- Cousins frequentist approach 22|| to determine the statistical sig- 
nificance of our measurement. In order to form confidence intervals, we use the A^ and 
S n7T distributions of the results of fits to MC pseudo-experiments for various input values of 
Ann- and S^ir- The distributions incorporate possible biases at the boundary of the physical 
region as well as a correlation between A wn and S n7T ; these effects are taken into account 
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by this method. The distributions are also smeared with Gaussian functions that account 
for systematic errors. The details of the method used to obtain the confidence intervals are 
described in Appendix [X] Figure |H1 shows the resulting two-dimensional confidence regions 
in the A^ vs. S nn plane. The case that CP symmetry is conserved, A n7T = S wn = 0, is ruled 
out at the 99.93% confidence level (C.L.), equivalent to 3.4<r significance for Gaussian errors. 
The minimum confidence level for A wn = 0, the case of no direct CP violation, occurs at 
(S n7T , A n7T ) = (—1.0,0.0) and is 97.3%, which corresponds to 2.2a significance. 

If the source of CP violation is only due to B — B mixing or AB = 2 transitions as 
in so-called superweak scenarios [111 EH' then (S nn ,A n7r ) = (— sin20i,O). The C.L. at this 
point is 98.1%, equivalent to 2.3a significance. 
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FIG. 8: Confidence regions for A^n and 



IX. DISCUSSION 

Using the standard definitions of weak phases 0i, 2 , and 3 , the decay amplitudes for 
B° and B to 7r + 7r _ are 

A(B° -> 7T+7T-) = -(\T\e i5T e^ 3 + \P\e i5p ), 

A(B° -> 7T+7T-) = -(|T|e i5T e"^ 3 + \P\e l5p ), (4) 
where T and P are the amplitudes for the tree and penguin graphs and 6t and Sp are their 



strong phases. Here we adopt the notation of Ref. 25] and use the convention in which 
the top-quark contributions are integrated out in the short- distance effective Hamiltonian. 
In addition, the unitarity relation V* b V u d + V* h V c d = —V£Vtd is applied. Using the above 



18 



expressions and 02 = 7r — 0i — 03, we determine 



_ 2 ^ l + lP/T|eW3) 
™ l + |P/T|e^3) • W 

Explicit expressions for S nn and are 

^tttt = [sin20 2 + 2|P/T|sin(0i — 2 )cos(!> 

A n7T = -[2|P/T|sin(0 2 + 1 )sin5]/ft, 
ft = 1 - 2|P/T|cos5cos(0 2 + 0i) + \P/T\ 2 , (6) 

where 5 = 5 P - S T . We take -180° < 5 < 180°. When A W7T is positive and 0° < fa + 2 < 
180°, 5 is negative. 

Recent theoretical estimates prefer \P/T\ ~ 0.3 with large uncertainties 2^, 27. liii 129^. 
Figures Efa)-(e) show the regions for 2 and S corresponding to the 68.3% C.L., 95.5% C.L. 
and 99.73% C.L. region of A n7T and S n7T (shown in Fig. |HJ) for representative values of \P/T\ 
and 0i. Note that a value of (S n7r ,A nn ) inside the 68.3% C.L. contour requires a value of 
\P/T\ greater than ~0.3. 

The allowed region is not very sensitive to variations of 0i within the errors of the 
measurements, as can be seen by comparing Figs. HJa), (c) and (e). The range of 02 that 
corresponds to the 95.5% C.L. region of A^ and in Fig. [H] is 

78° < 02 < 152°, 

for 0i = 23.5° and 0.15 < \P/T\ < 0.45. The result is in agreement with constraints on the 
unitarity triangle from other measurements |3l| . 



X. CONCLUSION 

In summary, we have performed an improved measurement of CP violation parameters 
in B° — > 7r + 7r~ decays. An unbinned maximum likelihood fit to 760 B° — > tt + tt~ candidates, 
which contain 163t23(stat) it + it~ signal events, yields A n7T = +0.77 ±0.27(stat) ±0.08(syst), 
and S nn = —1.23 ±0.41(stat) lam^syst), where the statistical uncertainties are determined 
from MC pseudo-experiments. This result is consistent with our previous measurement 0| 
and supersedes it. We obtain confidence intervals for CP-violating asymmetry parameters 
A n7T and S n7r based on the Feldman- Cousins approach where we use MC pseudo-experiments 
to determine acceptance regions. We rule out the CP- conserving case, A^ = S nn = 0, at 
the 99.93% confidence level. 

The result for S n7T indicates that mixing- induced CP violation is large, and the large 
A n7T term is an indication of direct CP violation in B meson decay. Constraints within the 
Standard Model on the CKM angle 02 and the hadronic phase difference between the tree 
(T) and penguin (P) amplitudes are obtained for \P/T\ values that are favored theoretically. 
We find an allowed region of 02 that is consistent with constraints on the unitarity triangle 
from other measurements. 
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FIG. 9: The regions for 02 and 5 corresponding to the 68.3%, 95.5%, and 99.73% CL. regions of 
A™ and S nn in Fig.EJfor (a) fa = 25.9°, |P/T|=0.3, (b) fa = 23.5°, |P/T|=0.15, (c) fa = 23.5°, 
\P/T\=0.3, (d) fa = 23.5°, |P/T|=0.45, and (e) fa = 21.3°, \P/T\=0.3. The horizontal dashed 
lines correspond to fa = 180° — fa. 
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APPENDIX A: MC PSEUDO-EXPERIMENTS AND CONFIDENCE REGIONS 

We use ensembles of Monte Carlo (MC) pseudo-experiments to determine the signifi- 
cance of our measurement and obtain confidence regions. They are also used for various 
crosschecks. Each pseudo-experiment consists of events that are generated with the nomi- 
nal PDFs, which are incorporated in Eq. EJ Since the parameters in the PDFs are derived 
from large-statistics control samples and sideband events, the pseudo-experiments precisely 
reproduce At distributions that are consistent with data. In particular, they are free from 
possible discrepancies between data and GEANT-based detector simulation. 

To generate each event in a pseudo-experiment, we first choose one LR-r region m ran- 
domly from a population that is based on the regional event fractions obtained from data. 
We then generate AE and Mb c values with distributions that are determined by the event 
fractions g™ n , and g^, which are listed in Table ITT1 The values of the probability func- 
tions f™, fx n and (in Eq. EJ) are determined from the AE and M bc values. We randomly 
choose an event type, nn, Ktt, qq or outlier, from a population based on f^, f^ n , f'^, and 
the outlier fraction f \. We generate q, At and resolution parameters according to the PDF 
of the selected event type. 

We repeat this procedure until the number of events reaches the observed number of 
events (760 events), and perform an unbinned maximum likelihood fit to obtain x^ nn and 
%Smr, which are the fit results and should be distinguished from the true (input) values A W7r 
and S nn . To account for the systematic error, each fit result is further modified by an amount 
determined from a Gaussian variation. We test the entire procedure using GEANT simula- 
tion, and find that distributions of x^ n7T and Xs-k-k obtained from the GEANT experiments 
are in good agreement with those from pseudo-experiments, when the resolution functions 
in the PDF are extracted from a lifetime fit to the GEANT data. We also verify that there 
is no fit bias as shown in Fig. 
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FIG. 10: Mean values of fit results vs. input values of MC pseudo-experiments for (a) A n7T and 
(b) iSj,-,,-. The solid lines are linear fit results. 
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We adopt the Feldman- Cousins frequentist approach |22j, which is based on the likelihood- 
ratio ordering principle, to obtain the confidence regions that are shown in Fig. |H] 32]. In 
the following, we first illusrate how we can obtain 1- dimensional confidence intervals for A^ 
with S n7T set to zero; intervals for are obtained in a very similar way. We then explain 
the method used for the determination of the two-dimensional confidence regions for A n7T 
and S n7r , which is an extension of that for the 1-dimensional case. 

We generate 10,000 experiments for 317 sets of (A™, S n7T ) values that cover the entire 
physical region. The fit to each set of experiments yields an xj^^ distribution that depends 
on the input A n7T value. To account for this dependence, we use a PDF for x An7T that consists 
of two Gaussian functions whose parameters depend on A^: 

P(x An7T \A n7T ) = f A ■ G(x Anw ] mi, cri) 
+ (1 - f A ) ■ G(x A7Tn ; m 2 , a 2 ), 

where G(x; m, a) represents a Gaussian function with mean m and standard deviation a, 
and f A , mi( 2 ), and o"i( 2 ) are polynomials of A^. The explicit expressions for f A , mi, m 2 , <J\ 
and a 2 are 

fA = 0,1 + a 2-4 7T7T 5 

mi = a 3 + a^A^, 
(J\ = «5 + o,qA W7T 2 , 

a 2 = an + a 12 A 

7T7T5 

where the 12 free parameters (a iy % = 1, 12) are determined from an unbinned maximum- 
likelihood fit to the x AlT7T distributions. Figures PHT a) and (b) show the distributions and 
the x An7T PDF for the cases (A n7r ,S nn ) = (0,0) and (Ann, St™) = (1)0), respectively. The 
PDFs are in good agreement with the distributions of pseudo-experiments in both cases. 
The acceptance region [x A i, x A2 ] for a given A^ and a confidence level a is defined by: 



a 



f x A2 

I dx P (jx ^n \Aji 



We adopt the likelihood-ratio ordering principle to determine x A \ and x A2 . Using the 
likelihood-ratio 

LR(x An7T \A 7Tn *) = P{x Al m | A-KTr) /P(x An7T | vAbest ) ) 

where Abest gives the maximum P value for a given x AlT7T , we require 

LR(x Ann \A n7T ) > Li?(^i|^4 7r7r ) = LR(x A2 \ A n7T ) 

for any x AlTn in [x^i,x^ 2 ]- Figure H21 shows the resulting confidence belts for A n7T . For a 
given measurement x An7T , a confidence interval at a confidence level a is obtained from the 
figure. 

The procedure to obtain the 2-dimensional confidence regions for A n7T and (Fig. |SJ) 
is an extension of the method described above. We use the following PDF: 

P\ X AjTn , XSnn {Aim , Sim ) 

nm m A\i Al) ' GySmr] m Sli &Sl) 

+ (1 - fAs) ■ G(A Kn ] m A2 , a A2 ) ■ G(S n7r ; m S2 , cr<? 2 ), 
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FIG. 12: Confidence belts in the A nn versus xa-k-k plane for a = 0.683, 0.955, and 0.9973 in the 
one-dimensional case. The dashed line corresponds to A^-n = xa-k-k- 

where /as, m Ai(2), m si(2), &ai(2) an d < J si(2) depend both on A n7r and S n7r . There are 27 
free parameters that are determined from an unbinned maximum-likelihood fit to the (x^ n7T , 
xs-k-k) distributions. We find that the PDFs represent the distributions of x_Akk and xs-kk 
very well for the input A wn and S n7T values over the entire physical region. An acceptance 
region Q at a confidence level a is also defined in a similar way to that for the 1-dimensional 
case: 
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where likelihood-ratio ordering is used. Using the requirement 

LRiyXj^-n j ■KSim \A nn > ^ittt ) ^ -£/-R( + 0.77, 1 .23 1 A nn , S nn ^j 

which corresponds to an acceptance region with our measurement (xAiw> x Smt) = 
(+0.77, —1.23) at its boundary, we scan the physical region in the A^^-S^ plane and cal- 
culate a confidence level a for each input point (A n7T ,S n7T ) to obtain the confidence regions 
shown in Fig. |HJ 

APPENDIX B: SOURCE OF SMALL MINOS ERRORS 

The Feldman- Cousins approach with acceptance regions determined from MC pseudo- 
experiments, which is described in Appendix |XJ is applicable to a wide range of analyses. 
On the other hand, care is needed when using experimental MINOS errors for the confidence 
interval calculation, as mentioned in Sec. IVIH Al In particular, difficulties may arise when 
the number of events is not large and the true values of physical parameters are located close 
to a physical boundary. In such small number of events can have a large influence 

on both the size of the MINOS errors and the shape of the log-likelihood ratio curve. The 
likelihood function for some events may become negative when the fit parameters are beyond 
the physical boundary. 

The observed features of the MINOS errors arise when there is an event that restricts the 
fit parameters in or close to the physical region, while the fit to all the other events gives 
a maximum likelihood that is located outside the physical region and is not allowed by the 
aforementioned restrictive event. For example, in this fit the removal of such a restrictive 
event results in an S n7T value that is more negative than S n7T = —1.23 (further from the 
physical boundary). In this case, the log-likelihood ratio curve is deformed by inclusion of 
the restrictive event, even if the curve before the inclusion is well-described by a parabola. 
The sizes of the MINOS errors also become small. 

We investigate this type of single-event fluctuation and its relation to the size of the 
MINOS errors with MC pseudo-experiments. For each experiment, we repeat the fit by 
removing each event in turn. The event that creates the largest difference in is tagged as 
the restrictive event and the change produced by the removal of the restrictive event, AS^, 
is recorded. When we choose the point of maximum likelihood at the physical boundary 
(A n7T , S W7r ) = (+0.57, —0.82) as the input for the MC pseudo-experiments, we obtain the 
average values of A«S T7r as a function of the positive error of «S 7r7r shown in Fig. O The 
correlation between the size of the error and the single-event fluctuation is evident. 

In our data, we have one event that has a large effect on the sizes of MINOS errors. 
The removal of this event from the fit gives S niT = — 1.9111q1 3 an d A^ = 0. 64^20, where 
the errors are MINOS errors; S nn is shifted to a more negative value (AS n7T = —0.67) and 
the MINOS error increases. This event has qr = —0.92 which is close to unambiguous B 
flavor assignment and corresponds to a very small wrong-tag probability. In addition, this 
event has AE = —0.01 GeV, and LR = 0.98, which corresponds to small B° — > K + tc~ and 
qq background probabilities. For this event, At = —3.8 ps and, thus, sin(AmrfAt)^ — 1. 
According to Eq. EH this event has a negative likelihood value at negative S nn values beyond 
~ —1.5, where it truncates the log-likelihood ratio curve. As a result, the negative MINOS 
error for the entire event sample is restricted by this single event. 

As shown in Fig. [T31 the observed single-event fluctuation AS n7T = —0.67 is consistent 
with the expectation from the MC pseudo-experiments if the positive error of is ~ +0.24, 
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FIG. 13: Single-event fluctuation versus the positive MINOS error on S n7T . The dashed lines 
indicate the observed AS^n and S wn positive MINOS error. 



which is the case for our data. A similar study for input values of S n7T and that are well 
within the physically allowed region indicates that this behavior occurs much less often. 
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